setwd("D:/作业/R语言")
load("D:/作业/R语言/class5_volcano.RData")
a=prostat
a$group=ifelse(a$FC>1.2&a$P<0.05,"Up",
               ifelse(a$FC<1/1.2&a$P<0.05,"Down","Not sig"))
ggplot(a,aes(x=FC,y=-log10(P)))+
               geom_point(aes(color=group))

  